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Abstract 


9mn 


one-dimensional Green's function 


A previously reported closed form solution is 
expanded to aetermine effects of isentropic mean 
flow on mode propagation in a slowly converging- 
diverging duct - a circular cosh ouct. On the 
assumption of uniform steady fluid density, the 
mean flow increases the power transmission coeffi- 
cient. The increase is directly related to the 
increase of the cutoff ratio at the ouct throat. 
With the negligible transverse gradients of the 
steady fluid variables, the conversion from one 
mooe to another is negligible, and the power 
transmission coefficient remains unchangeo with 
the mean flow oirection reversed. With a proper 
choice of frequency parameter, many different 
modes can be made subject to a single value of the 
power transmission loss. The paper also describes 
a systematic method to induce the effects of the 
gradients of the steady fluio variables. 


h 

I 

'^m 

K 

k 

k_ or k,(,p) 


k+ ,k,(,n) 
Lp»L^,Lx»L5 


radial variable in cylindrical 
coordinate system 

acoustic power flux 

Bessel function of order m 

total acoustic power across duct 

free space wave number w/Cq 

propagation constant of a mode or 
(m,n) mode in inlet 

propagation constant of a mode or 
(m,n) mode in exit 

linerar operations 


Nomenclature 

A(x),A(x) functions, Eqs. (61), (83) 


a 


Bm(x).Bo(x) 


b(x) 


effective length of converging- 
diverging section 

functions determining duct shape 
with or without mean flow, 

Eqs. (45), (48) 

duct radius varying along duct 




auct radius at throat (x = 0) 


b_ 


radius of left hand side uniform 
duct (inlet) 


M{x) 


Mo 

M_ 

M+ 

m, m' 

n, n' 

A 

n 

p 


Mach number of mean flow of con- 
stant steady density fluid 

Mach number at duct throat 

Mach number of uniform flow in 
i nlet 

Mach number of uniform flow in exit 
circumferential mode numbers 
radial mode numbers 
unit vector normal to surface 
steady pressure in flow 


b+ 


radius of right hand side uniform P^ 
duct (exit) 


steady pressure in flow of constant 
density 


Co 

Di , D2 , D< , 
D> ♦ Dp , D-t 


F 

G 


constant sound speed in isochoric 
flow 

sound speed in still fluid 
constants 


Pa 

P 

Qmn 


hypergeometric function 
Green's function 


q 

R 


deviation P - Pc 

acoustic pressure 

normalized Bessel function corre- 
sponding to (m,n) mode, Eq. (72) 

parameter, Eq. (56) 

acoustic power reflection 
coefficient 


■partially based on Consulting Reports to 
Pratt & Whitney Aircraft (Oct. 30, 1974 
and May 20, 1975). 

♦Aerospace Engineer; Member AIAA. 
♦♦Professor of Physics and of Aeronautics 
and Astronautics. 


R(x) or R(r) axial wave function, Eqs. (30), 
(40), (41) 

r radial variable in local spherical 

coordinate system (dr « dx) 

SCx,t) interaction term, Eq. (17) 


S^(x) temporal Fourier transform of ^ 

S(x,t) 

p 

T acoustic power transmission 

coefficient Po 


see Eq. (12) 

acoustic perturbation of density 
spatially varying steady density 


t 


time 


Pc constant steady density 


TL 


U 

Uc 


Ua 

IT 

V 

w 


X's 

X 


x,Xq 


z 


a or a^n 
6(x) 

r 


acoustic power transmission loss, 
Eq. (118) 

steady velocity of isentropic flow 

steady velocity of constant density 
fluio flow 

deviation IT - 

acoustic velocity field 

a constant, Eq. (57) 

a parameter, Eqs. (^5) and (46) 

w with M = 0 

see Eqs. (105) to (109) 

axial coordinate 

three dimensional coordinates 

dimensionless axial coordinate, 

Eq. (55) 

eigenvalue corresponoing to a ouct 
mode or (m,n) mooe 

b(x)/bo 

b±/bo 

Gamma function 


P5 deviation po - pc 

a a parameter, Eq. (66) 

T * b+/b_ 

♦("x) a Fourier transform of acoustic 

velocity potential 

4> azimuthal angle variable 

X see Eq. (13) 

’i'(x,t) acoustic velocity potential 

axial wave function, Eq. (41) 
u) angular frequency 

Introduction 

In a previously reported investigation of 
mode propagation in a nonuniform duct,l a closed 
form solution was obtained for a particular class 
of converging-diverging ducts - circular cosh 
ducts . These ducts are composed of two asymp- 
totically uniform ducts which are smoothly coupled 
through a converging-diverging section. The two 
uniform duct elements can differ from each other 
in cross sectional area. The ouct shape can be 
adjusted by means of three ouct parameters to 
produce a variety of circular ducts of practical 
interest. The analysis, which was previously 
developed for the case of no mean flow, is here 
expanded to include mean flow effects. 


Y specific heat ratio 

Yq mode cut-off ratio at ouct throat, 

Eq. (120) 

Y mode cut-off ratio at inlet, 

Eq. (119) 

6(x-Xq) Dirac delta function 

Kronecher's delta 

c a parameter, 1 for h = 0, Eq. (45) 

c a frequency parameter, Eqs. (53) 

and (122) 

j) normalized coordinate variable, 

Eq. (32) 

0 polar angle in local spherical 

coordinate system 

00 . half cone angle of duct segment. 

Fig. 1 


The mean flow under consideration is sub- 
sonic, and is assumed to be isentropic. Viscous 
effects are excluded from discussion. Numerical 
results will be presented only for cases of low 
Mach numbers, for which the steady density and 
pressure hardly vary from one location to another 
in the duct. However a systematic method is pre- 
sented to treat finite deviations of the steady 
fluid variables from the flow of constant density. 

Wave Equation 

As customary, the fluid variables are first 
decomposed into steady and fluctuating parts. The 
steady part will be further decomposed to treat 
its spatial variation. The steady density, pres- 
sure and flow velocity are denoted respectively 
by pqi P 4nd U; and the fluctuations by p, 
p and u. Retaining only the first order terms 
of the fluctuations, one obtains from the con- 
tinuity equation 

(It ^ * ^)p yPv • u + ^ 7 p^ • u + (^ • 0)p = 0. 


V 


V 


a parameter, Eq. (47) 
frequency parameter, kbo/a 


( 1 ) 


And from the Euler's equation 


(it * iJ • v)u (u • v)(J + |- 


5p . ^ p = 0 . 
0 ^^0 




With the potential flow condition v x u, one sets 


( 2 ) 


Here y is the specific heat ratio of the fluid, 
ana we have used the adiabatic relation between 
density ana pressure 


yP 

p = ^p. 
^0 


Also used are the zeroth order relations 


^ • (PqO) = 0 , 




v)() + vP = 0 


(3) 


(4) 


(5) 


The steaay variables are, in general, func- 
tions of spatial coorainates in a nonuniform duct 
with mean flow. The functions may be compli- 
cated. However, in a subsonic flow, the major 
portion of the steady variables may assume simple 
coorainate dependence, ana the deviations from the 
simple dependence can be small. We now decompose 
the steady variables further as follows: 


^ 6 ’ 


P = Pc " P^- 


ff = Oc ^ 0 ^. 


( 6 ) 


(7) 


( 8 ) 


Here pc. Pc snd Uq are a new set of 
steady variables which are simple functions of 
spatial coordinates, and p^, P^ and O 5 
are the deviations. In the present analysis, 

Pc and Pc are chosen to be constant and set 
equal to the steady density and pressure in the 
left side uniform ouct element (x < 0 , or the 
inlet). The variable Oc is the velocity of 
the constant density flow, and therefore, it is 
subject to the isochoric flow condition 


(J = 0 
c 


(9) 


On inserting Eqs. ( 6 ) to ( 8 ) into Eqs. (1) 
ana ( 2 ), one obtains 


(It 




P + 


u + C = 0, 


( 10 ) 


(It 




+ i- Vp + 7x = 0, (11) 


where 


C = V 


(pO, 




u + u 


vP 


6 ’ 


( 12 ) 


u = vv(x,t). 


(14) 


where v is the acoustic velocity potential. 

With the substitution of this equation, Eqs. (10) 
and ( 11 ) are combined to yield the equation 


2 

V f - 


1 

7 


(|_+ Oc • v) f = S(J,t) 


(15) 


where 




( 16 ) 


Equation (15) is, strictly speaking, a homogeneous 
equation for the acoustic field. It is, however, 
convenient to treat it as an inhomogeneous wave 
equation with the interaction term S. With 
S = 0, Eq. (15) is the homogeneous wave equation 
in a mean flow of constant density and pressure. 
The term S represents the interaction between 
the acoustic field and the mean flow field due to 
the spatial variations p§, P 5 and O 5 of 
the latter. Since no acoustic energy is created 
or annihilated, the interaction accounts for the 
elastic scattering of sound by the spatial non- 
uniformity of the steady fluid variables. 

Equation (15) is linear for the acoustic 
field. Thus the temporal Fourier transform yields 


7^<>(x) 


1 

7 


iO • v)^ *(x) = S (x) 

C 0) 


(18) 


Here 


♦ (x) = 


,t) e’"’^ dt. 


^ r ..1 

SJJ) = ^ {(-i. - Dc . v)xJJ)^^U5)j. 


(19) 


( 20 ) 


with 


X (x) = — ^ / x(x, 

t (J) = ^ / t(x. 


t) dt. 


e dt . 


( 21 ) 


( 22 ) 
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Homogeneous Solutions 

We consider solutions to Eq. (18) first with 
the interaction term neglected. The wave equation 
is then 



( 28 ) 


1 d 
Sin e "?? 


sin 0 


dj 


+ 





(29) 


where the subscript H has been used to indicate 
the homogeneous solutions. The homogeneous solu- 
tions serve two purposes: First, with boundary 

conditions properly satisfied, they describe mode 
propagation in the duct with mean flow of uniform 
density, which can be adequately described by the 
variables p^. Pc Oc- It is likely 

that potential flow solutions will be close to 
these variables for low Mach numbers. Secondly, 
the homogeneous solutions can be used as basis 
functions to construct the Green's function. With 
the help of the latter, Eq. (18) can be system- 
atically solved by means of the iteration- 
perturbation method for cases when the deviations 
of Steady fluid variables are finite but small. 2 

Let us consider Eq. (23) in a circular duct 
which is specified in terms of the radius b(x), 

X being the axial coordinate variable. As in 
Ref. 1, we consider first a small duct segment and 
the local spherical coordinates (r,e,4>) as 
illustrated in Fig. 1. The azimuth 4 >, the 
angle variable around the auct axis, is not shown 
in the figure. The duct segment is so short that 
it may be regarded locally conical. Thus the 
slope, b = ob/dx, of the duct wall remains con- 
stant within the segment. The origin of the local 
spherical coordinate system is located at a point 
on the duct axis such that the coordinate surface 
e = 0g (cone) tangentially contacts the duct 
wall in the segment. The half cone angle 0 q 
is related to the wall slobe b' as 

b' = tan 0^. (24) 

With the spherical coordinates, the isochoric flow 
velocity is chosen as 



where m and a are the separation constants, and 



The solution to Eq. (28) is 


H = e”"’ (31) 

with m being an integer. As for Eq. (29), it is 
convenient to introduce the normalized coordinate 
variable n, defined as 


sin 0 


for ^ 0 


= ^ for e„ = 0. 

D 0 


(32) 


where h is the radial variable in the cylindri- 
cal coordinate system. The value of n varies 
from 0 to 1 in the duct. On the duct axis, n = 0; 
and on the duct wall, n = 1. Equation (30) is 
written in terms of n as 


c 


c 



(25) 


(> 


sin^0 






- 2 sin^0^n^) 


dn 


where « is a constant with the sign being 
opposite for converging and diverging sections, 
and r is the radial unit vector. Equation (23) 
can then be written for the spherical coordinates 
as 

ija 23. 1 3.. 

^ I ar 3r sin 0 30 ® 30 



For a duct of slowly varying cross section 
(|b'| « 1), sin 0Q = b'. Neglecting terms 
involving (b')^ and higher orders, one can write 
Eq. (33) as 



This equation can be separated into three ordinary 
differential eouations. With the substitution 


J(e) h(4)) (27) 



This is the Bessel equation of order m, and the 
physically acceptable solution is 

J 'v, JJan) (35) 


k 


one obtains, from Eq. (26) 


where Jm is the Bessel function of order m. 
From the boundary condition of the hard duct wall, 
one determines values of a as follows: 


a 


n = 0, 1, 2 . . 


= a 


mn’ 


(36) 


where a^n is the n-th zero of dJu,(x)/dx. 

Some remarks should be made on the coordinate 
variable n. Unlike the variable e, n is not 
a local coordinate variable, but it is valid 
throuoh the entire duct. To see this, let us con- 
sider two small duct segments 1 and 2, respective- 
ly at X = XI and X2* The half cone angles 
of the respective segments are denoted by 
00,1 ana 0o,2- Assume 0 q,i > 0q, 2- 
The variable 0 varies from 0 to 6q i in seg- 
ment 1, and from 0 to 0 q, 2 i>^ segment 2. A 

value of 0 greater than’ ©o,2 i^ss 
than 0 Q 1 , represents a coordinate surface in 
segment I, but not in segment 2. Furthermore, two 
coordinate surfaces represented by a value of 
0 less than 0o,2» respectively in the two 
duct segments are not, in general, equivalent to 
each other. On the other hand, n varies from 0 
to 1 for all the duct segments. The locations re- 
presented by a value of n are in one-to-one 
correspondence for two different duct segments. 

A coordinate surface is uniquely determined as the 
locus of n = ni, ni being a constant, and 
the coordinate surface spans all the duct seg- 
ments. A coordinate surface with a value of n 
does not intersect with any other coordinate 
surface with a different value of n. The duct 
wall is a coordinate surface corresponding to 
n = 1. We also note, from Eqs. (35) and (36), 
that an eigenfunction assumes the same value on a 
coordinate surface n = ni- A mode is in one- 
to-one correspondence with an eigenfunction. Thus, 
with the use of the variable n, a mode is defined 
in a nonuniform circular duct with isochoric flow 
if |b'| remains small (Fig. 2). 

Equation (30) governs the propagation of a 
mode along the duct! The equation is first trans- 
formed from the local coordinate r to the duct 
coordinate x. lo this end we use the relations 


sin 0^ = b' 

(37) 

b 

(38) 

II 

(39) 


Using Eqs. (37) to (39), one can write Eq. 
(30) as 


(1 - M^) (1 + M^) + ZIIcm] 

dx^ L J 


dx 




( 40 ) 


Here we have used 


M(x) 




where bo and Mq are the duct radius and 
the Mach number at the duct throat. For a given 
b(x), M(x) is determined once Mq is known. 
Thus, for a given duct shape, Mq is sufficient 
to specify M(x). 


With the substitution 


R(x) 


*00 

b(x) Vl - M^(x) 



M(x) 

1 - M^(x) 



(41) 


Equation (40) is transformed into 



where terms like b'2 and b" have been neglected. 


Solutions in a Circular Cosh Duct 

Equation (42) was solved exactly in a closed 
form with a circular cosh duct when M « 0, in 
Ref. 1. With M ^ 0, it is convenient to introduce 
a function Bm(x) as follows: 


These relations are valid for all the duct seg- 
ments and the error is of order (b')^. In order 
to avoid misunderstandings, a remark should be 
made on Eq. (39): Eq. (30) governs the wave along 

a line that is the intersection of coordinate sur- 
faces 4» = constant and n = constant (not 

0 = constant). On this line, dr = dx vTTT^, 
and Eq. (39) follows. One should not attempt to 
derive a similar relation by regarding r as a 
function of h as well as x. Such a derivation 
might lead one to an incorrect conclusion, and the 
error would appear to be of order b' instead of 
(b')2. Note that the origin of the local coor- 
dinate system changes with x. It cannot be over- 
emphasized that R (as given in Eq. (30)) deals 
with the wave field variation along a line which 
is the intersection of coordinate surfaces 
4» = constant and n = constant. Finally, the 
use of the relations (37), (38), and (39) is very 
unlikely to lead one to any asymptotic paradox. 


Bm(x) 


(1 - M^)b2 


V 


2 M^(2 - M^) 
(1 - M^) 


(43) 


where v = kbQ/a. Note that this function de- 
pends on the Mach number and the frequency (v) as 
well as the duct shape (b(x)). With the use of 
Eq. (43), (42) is written as 




- B^(x) 

‘’o 


= 0 . 


(44) 


The discussions thus far apply to any nonuniform 
hardwall circular duct as long as b'2 and b" 
are negligible. 

To utilize the cosh duct solution, we set 


Bi^(x) = Bj^(o) + w |cosh^p sech^(ex/a - p) 

- sinh (2p) tanh (ex/a - p) - cosh (2p)| . (45) 

Here a is the effective length of the 
converging-diverging section, ana w, p ana e 
are parameters depending on M, v ana duct geo- 
metry. The parameters w and p are determined 
as 

w = {[B^(o) - B„(-)][B„(o) - B„(-)]}^/^ (46) 


P 



Bm(o) - 


( 47 ) 


Involved in these equations are M, v and b± 
the radius ratios of the duct. For fixed values 
of Mq 3^nd V, w and p are ^determined from 
the radius ratios e±, and vice versa. One 
notices from Eq. (47) that p is equal to zero 
for a symmetric duct (e+ = e-), regardless of 
values of Mq ana v. 


n I_c. 1 2 M^(-s)[2 - M^(-s)] 

^ [1 - M^(-s)]b^(-s) [1 - M^(-s)]2 

(51) 

and Eq. (45) is 

B|^.(-s) = B|^(o) + w jcosh^u sech^(es/a + y) 

+ sinh (2w) tanh (cs/a + y) - cosh (2y)| (52) 

The value of e is determined from Eqs. (50) to 
(52). There may exist more than one real solution 
for e; the correct one should be determined from 
physical considerations. A list of values of e 
are given in Table 1, for various values of Mq 
ana c, where c is a frequency parameter 
defined as 



The parameter e is not given in a closed 
form, but its value can be determined analyti- 
cally. For fixed values of Mq and v, the 
duct radius ratios e± have been used for the 
oetermination of w and p. The axial variation 
of the radius in the converging-diverging section 
will be used to determine e. If the value of e 
were kept constant, the duct shape would depend 
on M and v as can be readily seen from Eqs. 
(43) and (45). Thus, the value of e must be so 
determined that the duct shape may remain un- 
changed for different values of M and v. To 
this end, at a properly chosen axial location in 
the converging-diverging section, the ouct radius 
is required to be the same for all the allowed 
values of Mq and v. We choose Mq = 0 as 
the reference case. With Mq = 0, we set 
E = 1, and obtain from Eq. (43) 


As the table shows, e is close to unity in most 
cases of interest. Also notice that e takes the 
same value for all the modes at c = 0. The 
value of E at c = 0 are listed for various 
values of Mq and b_ in Table 2. 

We now return to Eq. (44). On inserting Eq. 
(45) into Eq. (44) we obtain 

+ j^q^ - v(cosh^p sech^z 

- sinh (2p) tanh z^J = 0 (54) 

where 


Bq(x) i B^(x) 


M 


= b“^(x), 

0 


( 48 ) 


and from Eq. (45) 

Bp(x) = 1 + |cosh^y^ sech^(x/a - 

- sinh (2yp) tanh (x/a - y^) - cosh (2yQ)j (49) 

where the subscript o indicates M = 0. Note 
that Bq, Wq and pq do not depend on v. 

The axial location is chosen preferably in the 
converging section, as x = -s. The duct radius 
at that location is determined from Eqs. (48) and 
(49): 


z = EX/a - y ( 55 ) 

= (ka/c)^ - (aa/cbp)^ ( 2 y)]. 

(56) 


V = (aa/ebjj)^ w ( 57 ) 

Equation (54) can be solved exactly. 2 The gen- 
eral solution is written as 

= D^4»^(z) + D24»2(z) (58) 

Here and D 2 are constants; ana 

and are linearly independent solutions. 


b(-s] 


■{ 


1 + 


2 

w^ cosh p^ 
0 0 


sech‘^(s/a + 






*1 = A(z)Flcp?2>43J 






( 59 ) 


Sinh ( 2 yp) tanh (s/a - y^) - cosh (By^) ( 50 ) 

With Mq ^ 0 , Eq, ( 43 ) is written for x = -s as 


-ik^a/E 

♦2 = A(z) (1 + e^") 


6 


*^(^2 ■ ^3 * ■ ^3 ^ ■ ^ 3 ’ )’ 


(60) 


where F is the hypergeometric function, and 
i (k^-k_)za/2c 


A(z) = e 


(2 cosh z) 


Co = 1 - ik+a/€ 


1/2 


:^+k_)a/2£ 

(61) 

eT '' 

(62) 

"). 


.i}/2 

(63) 

.))/? 

(64) 


(65) 


(66) 


0 = (4v cosh u - 1) 

The constants ana D 2 are to be 
determined from bounoary conditions at x = 

The asymptotic forms of the hypergeometric func- 
tion are used to find 

i ) at X = +CO 

-ik^pa/e ik^x 


= e e 


ik^pa/e -ik^x 
4»2 = e e , 


i i ) at x = -* 

r(c 3 )r(ci ^2 ^ pa/e+ik_x 

♦1 ® 

r(c 3 )r (;3 - - ^ 2 '^ ik_(ia/c-ik_x 


(67) 

(68) 


r(C3 - Ci)r(c3 - C 2 ) ® 


(69) 


r(2 - C 3 )r(c^ ^2 “ ^3^ -ik_pa/e+ik_x 

'*'2 ^ r ( c 2 - ^2 ^ 


r(2 - C 3 )r(c 3 ~ - ^ 2 ^ ik_pa/e-ik_x 

r(i - c^Tm - q) ® 

where r is the Gamma function. 


Mode Propagation in Circular Cosh Ducts 
With Low Mach Number How 

As mentioned earlier, the deviations p^, 

Pj and 0^ are negligible in cases of low 
Mach numbers (Mq < 0.5). In these cases, the 
solutions can be obtained as in Eq. (58), with the 
boundary conditions to be satisfied. Note that 
these solutions are the zeroth order solutions to 
Eq. (18) and can be used to obtain more accurate 
solutions. 

Let us begin with the incident wave that is 
composed of the (m,n) mode coming from x = 

In the incident side uniform duct (x << 0), 




im4> 1 


X exp 


b_ 

^i^-kM_/(l - M^) + (71) 


where 




1/2 


k^*^ = k\ 

mn * 


(72) 


(73) 


(74) 


mn 


Corresponoingly, we have 


= e 


(75) 


With mode conversion neglected, the solution ^ 
is given at x ±* in the form 

ik^^x — ik^^x 

xll!!- = e + Dj. e , (76) 


ikitL 


lim i/v#\ r r\ 

X^oo ® 


(77) 


The second term of Eq. (76) represents the re- 
flected wave, and Eq. (77) the transmitted wave. 
The constants in Eq. (58) are readily determined 
by comparing the asymptotic expressions in Eqs. 
(67) to (70) with Eqs. (76) and (77): 


= 0 


Oj.e 




r(c^)r(c2) 


(78) 


(79) 
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The coefficients and Dt are 


D = e 
r 


r{c^)r(c2)r(c3 - 

r(Cl ^2 " " ^2^’ 


0, = e 


r(C3)r(c^ + ” C 3 ) 


(80) 

(81) 


In Eqs. (78) to (81) the parameters w, e, ci» 
C2 and 53 contain the constant a: it 

is to be replaced by the eigenvalue amn« The 
solution is 

♦ = Q^^(n) e^"’* b-l(l - 

(f- 

A(z) r^CpC2,C3; (1 + j 


1 V 
= 27 A " 


n 

Here g^n is the one-dimensional Green's func- 
tion. On inserting Eq. (86) into Eq. (85) and 
using the orthogonality and the completeness of 
the eigenfunctions, one obtains 





(87) 


X exp 


dx . (82) 


where 


The homogeneous solution of this equation is none 
other than R given in Eq. (41). 

Consider the factor 4* first. For a cir- 
cular cosh duct, the homogeneous solutions 
and \|72 are given in Eqs. (59) and (60). 

The source term is non- zero only in the 
converging-diverging section. Thus, the Green's 
function is required to represent outgoing waves 
in the uniform duct elements. Consider a new set 
of homogeneous solutions 4»> ana 4<< which 
asymptotically behave as 


kM 


^ X + ik 

L 

♦^(x) l:; e . 

(88) 


(83) 

-ik(-) 



♦^(x) e . 

(89) 


x_ being a large negative value of x. Recall 
that n = h/b(x) and z = ex/a - y. 

Green's Function ana General Solutions 

When the interaction term is finite 
but small, Eq. (18) can be solved by means of the 
iteration-perturbation method. 3 To this end, 
the equation is first replaced by the integral 
equation 

♦ (J) = + y* G(X,J„)S^(J^) d\ (84) 

Here G is the Green's function to be discussed, 
and the function 4>(^) has been used for the 
homogeneous solution, to let the solution satisfy 
the boundary conditions. 

The equation for the Green's function is 


Comparing these equations with Eqs. (67) to (70), 
one can readily find 

= ♦!. (90) 

r(c3)r(c2 - ^3 l)r(c3 - S3 + 1) 


*< = *1 - 


r(sj)r(s2)r(2 - S3) 


(91) 


Corresponding to these functions, the homogeneous 
solutions to Eq. (87) are 

Rj_(x) = b“^ (l - M^) exp J ^ ^ dx^ ♦j.(x). 


^lii + iO^ • G(x,x^) = i(x - 


Jp) (85) 


(92) 


The solution can be obtained in terms of the homo- 
geneous solutions with proper boundary conditions 
as follows: 


!<<(x) = b ^ (l - M^) exp ^ik j dx^<i>^(x). 


(93) 
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The one-dimensional Green's function is then ob- 
tained as 


” X a(x 


R>(x) R^(Xq) for X > 


R (x) R (x^) for X < X 
< ' > o' 0 


s(0) = L 
0 ) S 

with the operator defined as 



( 102 ) 


(103) 


where a 1 s the Wronskian 


R r' 
< < 

R r' 
> > 


(94) 


(95) 


On inserting Eqs. (86) and (102) into Eq. 
(84), one obtains 

♦ (x) = »^°^(x) + — ^ 

"^2 IT 


s 

m ' n ' 




(104) 


The source as given in Eq. (20) is a 
linear function of acoustic field variables. The 
latter are not known until the t>roblem is solved. 
However, when the deviations, p^, and 
U5 are small, the source is small. Thus, the 
zeroth order solution «>(0) can be used for 
the acoustic field variables to compute 
The acoustic velocity is u(^) = V4>v0)^ and 
the acoustic pressure is obtained from Eqs. (10) 
and (11) with ^ and x neglected. 


where the incident wave is the (m,n) mode (see 
Eq. (82)), and 

^mn;m'n'^^) = 


1 /’2’' /*1 r -im'4'o 

’‘IT ‘^*0 J ‘*’‘ 0 ® Om'n'<\ 


P 


( 0 ) 



(96) 


9, 


m'n' 


(x.Xo)Ls e 


im(|> 




(105) 


where Lp is a linear operator defined as 

'■p ■ ' |(tl) 


For an axisymmetric flow, is independ- 
ent of 4>o thus one can readily obtain 


y - X y'" 

mn;m'n' ®mm' nn' 


(106) 


(97) 


nn ' 


= f \ Qmn'(''o^9mn'(’‘>’‘o) 


It fol lows that 


Ls Qa,n(''o)'^ 


( 0 ), 


Xo)- 


(107) 




(98) 




(99) 


Here and are linear operators 
containing the deviations of the mean fluid 
variables: 


Furthermore, if p^, and U5 are in- 
dependent of n and if n • (J^ = 0, one has 
then 


C- =^nn- 


= / dx g (x,x )L_ R 
nn y 0 ^mn' ’ 0' s 


( 0 ) 


(Xo). 


(108) 

(109) 


= (? . 0, - 0, 


^)Lp 


p I w ^b * ^ ^ 


VP. 


V, 



It follows that mode conversion, from one mode to 
another, occurs only if the transversal gradients 
of the steady fluid variables are nonnegl igible.^ 

(100) Power Reflection and Transmission Coefficients 

The acoustic power intensity in the isen- 
tropic flow is given by^ 

^ 2 2 \ 

T = <pu> + uf<p^>/p^c^ + 0 • <pu>/cj 


The source term is then 


+ pQ<(0 • u)u>, (lie) 
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where < > stands for the time average, 

total power across the duct is obtained as 

The 

K =y* T • n dA 

(111) 

where oA is the surface area element and n the 
unit vector normal to the surface. In the uniform 
duct element containing the axisymmetric mean 
flow, we have 

K = y I^h dh 

(112) 

with the axial component I^ of the power 
intensity given by 





(113) 

Let Ki, Kp and Kt denote the total 
powers respectively of the incident, the reflected 
and the transmitted wave. The power reflection 
and reflection coefficients, R and T, are de- 
fined as 

R = K^/K., 

(114) 

T = K^/K.. 

(115) 

Using Eqs. (71), (81), (76) and (77), one can 
obtain 

cosh [Tr(k_^ - k )a/e] + cosh (tio) 
^ ~ cosh lv(k^ _ k_^)a/ej + cosh (tio) 

(116) 

2 sinh (irk^a/e) • sinh (irk a/e) 

^ - cosh LTr(k^ + k^)a/ej + cosh {vo} 

(117) 


For the present calculation, Eq. (117) is used 
for T. The TL is plotteo as a function of the 
throat Mach number Mq or of the frequency 
parameter y- c. The parameter y- is 
the mode cut-off ratio referenced to the left hand 
side uniform duct element, and is referred to as 
the inlet cut-off ratio: 


Y_ = kb - (119) 

The inlet cut-off ratio is a convenient parameter 
character iz ing the modes contained in the incident 
wave. The sound generated in a fan duct often 
comprises many different modes. In such case, the 
mode distribution can be obtained as a function of 
the cut-off ratio. Another cut-off ratio useful 
for discussion is that referenced to the duct 
throat, given by 


If Yo > The mode is cut-on through the en- 
tire duct. On the other hand, if yq < The 
mode that is cut-on initially (y_ > 1), is cut- 
off in the converging-diverging section. Except 
for low frequencies (kbo < 3), the TL is ap- 
proximately 3 dB for Yo = independent of the 
duct geometry or the Mach number. The two cut-off 
ratios are related to each other as 

v_ = V-[(l-"o) (l-"o (121) 

The parameter c is defined in Eq. (53), and 
can be written as 



Note that R + T = 1. That is, the acoustic 
energy is conserved as it shoulo be in a linear 
analysis of the acoustic field. The terms neg- 
lected in the approximation woulo not create or 
annihilate the acoustic energy. They, if in- 
cluoea, might transfer the energy between mooes or 
between the reflected and the transmitted waves. 
Also note that Eqs. (116) ana (117) are similar to 
Eqs. (25) and (26) in Ref. 1. The mean flow 
oepenoence is incluoed only through the parame- 
ters k±, € ana a. The latter are even 
functions of M. Thus the power reflection and 
transmission coefficients are inaepenaent of the 
mean flow direction, the positive or the negative 
x-direction. 

Numerical results will be discussed in ternis 
of the acoustic power transmission loss (TL) in 
the circular cosh ducts. The mean flow is isen- 
tropic flow ana has a uniform steaay fluid den- 
sity. As mentioned earlier, a uniform steaay 
density is a good approximation for low Mach num- 
ber flow in inlet ducts. Numerical calculation 
includes the Mach numbers up to the value of 0.5 
at the duct throat. 

The TL is defined as 

' TL = -10 log^Q(T) (118) 


This parameter collapses the TL curves for many 
different modes. 

In Fig. 3, the TL oivided by kb_ is 
plotted as a function of y~ Tor the various 
values of Mq between 0 and 0.5, ana for 
3_ = 1.1 or 1.3. .The other duct parameters are 
a/b_ = 1, c = 1. On the top of the figure are 
shown the arrow marks, each pointing to a value 
of Y- which corresponds to yq = i w*Th the 
given value of Mq (cf. Eq. (121). The (8,5) 
mode has been used for the calculation. However, 
except in the vicinity of Yg = The results 
may be used for modes with eigenvalue from 8 to 
25, within an error of 1 dB. As expected, the 
TL decreases with the increasing value of y- 
(see Ref, 1). The decreasing rate of the TL is 
faster for larger Mach number. For a fixed value 
of Y-» The TL is smaller for larger Mach num- 
ber. Since, TL/kb_, is plotted in Fig. 4, TL 
itself is proportional to kb_ for a fixed 
value of Y-- IT follows that, for a mode dis- 
tribution given as a function of the inlet cut-off 
ratio, the converging-diverging auct contour re- 
duces the transmitted sound more effectively for 
cases when the distribution comprises modes of the 
larger eigenvalues. 
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In Fig. 4, the TL/kb_ is plotted as a 
function of Mq, for fixed values of y_. As 
Mq increases, the TL decreases very slowly in 
the beginning (Mq 2^ 0), and then decreases at a 
gradually increasing rate until Mq reaches a 
value for which yo = 1- mean flow depend- 

ence of the TL is directly related to the change 
of the cut-off ratio at the throat. As one can 
notice from Eq. (112), yq increases with in- 
creasing Mq for Y- held constant. 

In Fig. 5, the TL is plotted as a function 
of c. Each curve, for a value of Mq, in- 
cludes many modes with eigenvalues greater than 
5. The TL increases with increasing c, and 
becomes approximately 3 dB at c = 0 independ- 
ent of the Mach number and the duct parameters. 

The advantage of this presentation is that almost 
all the modes are subject to the same TL for a 
given value of c. The mean f^ow dependence of 
the TL in this figure appears misleading. How- 
ever, note that, with Mq increased, yo in- 
creases and kbo should be decreased to keep 
c unchanged (see Eq. (122)). 

Concluding Remarks 

In an attempt to improve the understanding of 
the acoustic characteristics of a fan duct system, 
mode propagation has been investigated in a par- 
ticular class of converging-diverging circular 
ducts - circular cosh ducts, with isentropic and 
inviscid mean flow. The ouct shape can be ad- 
justed by means of three duct parameters, and 
covers a wide range of coverging-di verging ducts 
of practical interest. 

On the assumption that the ouct cross section 
area varies slowly, an approximate wave equation 
has been derived. The equation is divided into 
two parts: (1) the homogeneous wave equation 

involving the mean flow of uniform steady density 
fluid, and (2) the interaction between the acous- 
tic field and the deviation of the steady fluid 
variables from the constant steady density fluid 
flow. The homogeneous equation is solved in a 
closed form, and the interaction term is treated 
by means of an iteration-perturbation method. The 
solution of the homogeneous equation is regarded a 
good approximation in cases of low Mach number 
mean flow. 
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TABLE 1. - VALUE OF e 



Mo 

Eigenvalue 





-0.4 

-0.2 

0 

0.2 


0.2 

4.20 


1.0020 

1.0014 

1.0007 



8.54 

1 . 0020 

1.0017 

1 

1.0010 



16.53 

1.0017 

1.0015 

1 

1.0012 

1.1 


23.27 

1.0016 

1.0015 


1.0012 

0.4 

4.20 



0.9920 

0.9845 




8.54 

0.9980 

0.9952 

1 

.9886 



16.53 

.9953 

.9937 

1 

.9903 



23.27 

.9944 

.9932 

t 

.9908 


0.2 

4.20 

1.0039 

1.0026 

1.0012 

0.9997 



8.54 

1.0025 

1.0019 

1 

1.0004 



16.53 

1.0019 

1.0015 


1.0008 

1.3 


23.27 

1.0017 

1.0014 

i 

1.0009 

0.4 

4.20 

1.0041 

0.9948 

0.9839 

0.9710 




8.54 

.9946 

.9895 

1 

.9778 



16.53 

.9897 

.9868 


.9808 



23.27 

.9880 

.9860 

+ 

.9817 


With the interaction neglected, a mode is 
preserved, and the acoustic power transmission 
coefficient increases with the mean flow in- 
crease. This mean flow effect is directly related 
to the change of the cut-off ratios; With the 
increasing mean flow, the cut-off ratio at the 
duct throat increases faster than at the uniform 
duct sections (inlet and exit). For given duct 
geometry and mean flow, the TL of many different 
modes can be made to collapse onto a single curve 
with a proper choice of frequency parameter (c). 
The mean flow effects remain unchanged with the 
flow reversed. The interaction term needs numeri- 
cal computations for its detailed discussions. 
However, the formal solution shows that a mode may 
be converted to others only if the transverse 
gradients of the steady fluid variables are not 
negligible. 


TABLE 2. - VALUE OF e AT c = 0 


Mo 



1.1 

1.2 

1.3 

1.4 

1.5 

0.1 

1.0004 

1.0006 

1.0005 

1.0004 

1.0002 

.2 

1.0014 

1.0015 

1.0012 

1.0005 

.9998 

.3 

1.0008 

1.0001 

.9987 

.9971 

.9954 

.4 

.9920 

.9875 

.9837 

.9806 

.9776 

.5 

.9438 

.9316 

.9271 

.9244 

.9223 
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DUCT SEGMENT 







’*M(r,0. 0))] 
M 


Figure 1. - Duct segment and local spherical coordinate system. 



Figure 2. - Coordinate surfaces (17 «= 0.345, 0.719, 1) and trans- 
versal shape of ( 1 , 1 ) mode at various axial locations. 



Figure 3. - Power transmission loss (TL)of (8.5) mode, divided by kb., 
versus Inlet cutoff ratio (y.). 


TL{dB)/kb. 



Figure 4. - Power transmission loss (TL)of (8.5) 
mode, divided by kb_, versus Mach number at 
throat (Mq). 




Figure 5. - Power transmission loss (TL) versus 
Each curve representing modes with eigen- 
value greater than 3. 
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